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We present an efficient algoritlim for calculating multiloop Feynman integrals perturbatively. 

1. Recently, a new method has been proposed to calculate Feynman integrals of multiloop diagrams perturbatively 
0] . Together with the solution procedure for graphical recursion relations developed in Ref . , this should ultimately 
lead to the completely automatized computer generation of perturbation expansions of field theories up to high orders. 
Such expansions are needed in all strongly coupled fluctuating field systems, for example those describing the critical 
phenomena close to second-order phase transitions So far, expansions have been limited to seven loops only 
^\ [ which are barely sufficient to yield critical exponents Q with an accuracy comparable to experimental data 
0^ ' In this note, we would like to show how the expansions proposed in |Q] can be performed most efficiently, such that 
0^ ■ they can be carried out on a computer to high orders in a limited computer time. 
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2. A basic Feynman integral with L loops, n internal lines, and E external momenta ki, . . . ,k£; has the form 



, with some powers a^, where are the momenta carried by the lines, and the integrations run over all loop momenta 
^ ' Pi . The line momenta are linear combinations of the loop momenta and the external momenta kj . 
[- — For simplicity, we have set all masses equal to unity. A Feynman integral with different non-zero masses can be 
reduced to ([|) by an appropriate rescaling of the line momenta qfc. Following Ref. [|], we view the integral (|l|) as a 
special case = (1) of the function 



to be calculated perturbatively via a Taylor series expansions in powers of k. 

It is the purpose of this note to point out that the simplest way to derive such an expansion is by rewriting each 
generalized propagator in the Schwinger parametric form 
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a(K-l)q^ 2 



(I + «q2)- r(a) 7o 
Then the integrals (p|) take the form 
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fc=l 

Collecting the L loop momenta pi and the E external momenta kj in single vector symbols 



p = (pi, . . . ,Pl), (5) 

fc=(kl,...,k£), (6) 



we rewrite 
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^[afe(l-K) + Kife]q2 = -p^Mp + p^M'k + -k^M"k (7) 



fc=i 

and complete the squares to 



Yyak{l-K) + Kife]q^ 2 ^ M-^M'kf M {p + M'^M'k) + -fc^ (m" - M'^M'^M') fc, (8) 

with symmetric matrices M and M". After a shift p ^ p — M^^M'k of integration variables, the pfe integrations 
become Gaussian, and we obtain 



(n \-LD/2 /-oo /.oo ~^k'^(M"-M''^M-^M')k 



where the matrices M, M', and M" depend on k and the through linear combinations 

CkiK,tk) = akil - k) + Ktk- (10) 

Although the entries of the matrix M depend on the routing of the loop momenta through the different lines, the 
determinant of M is invariant under changes of the routing, except for trivial relabelings of the ak- 
in order to derive the desired expansion of Iaf,{K) in powers of k, we expand the integrand on the right hand side 
of in powers of k, whose coefRcients are polynomials in the parameters ti, (i — 1, . . . , L). The t^-integrals can then 
all be performed using the formula 

dtt^e"* = r(7 + l). (11) 

For diagrams without external momenta, appearing in the perturbation expansions for the ground state of quantum 
field theories, (^ simphfies to 

O \~LD/2 poo nOO 

laM- ^ , / rfiiC-'e-'^-.. / dUtr'e-HdetM)-^/'. (12) 

T{ai)---T{ai) Jo Jo 

More general Feynman integrals than those in Eq. (|l|) may contain loop momenta pfc in the numerator of the 
integrand. These can be calculated with a simple extension of the above technique, by introducing "source terms" 
■ Pi into the exponents of (||) and (||), and appropriately differentiate the resulting K-expansion with respect to 
ji, which are set equal to zero at the end. 

3. As a first example, take the exactly solvable one-loop integral 

r,^ f d^p 1 ^ T{a - D/2) 

" J{2n)D{l + p2Y (4^)^/2r(a)- ^ 
Its K-generalized version can be expressed in terms of a confluent hypergeometric function, 
^a, , f d^p e"(''-i)p' 1 /D ^ D a{l - n)\ 
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with 



r(a) ' ' V 2 2 ' K J r{D/2) ^ ' V ' 2' k 



iFi(a;/3;z)^ V)^i-, (a), ^ = TT (a + r) (Pochhammer's symbol). (15) 
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In Ref. [Q this was calculated perturbatively via a Wick expansion. Here we use our general formula (|T^) for vacuum 
integrals. The number of loops is L = 1, and we identify 

qi — p, fli = a, ci = a(l — k) + Kt, M = 2(ci), detAf = 2ci. (16) 

Expanding (detM)~^/^ in powers of k, and performing the resulting integrals over t in Eq. (p^, we find directly the 
perturbation expansion for the loop integrals (|l3|) in any dimension D: 



(47ra)^/2 



D{2 + D)k^ D{2 + D){4: + D)k^ (2 + a) £> (2 + £>) (4 + Z?) (6 + £») 



8a 24a2 128a3 

(17) 



{6 + 5a)D{2 + D){i + D){6 + D)_(8 + D)K^ e 

h U[k 



960a4 

The expansion can easily extended any desired order. It agrees, of course, with what we would obtain from the exact 
expression (n4h via a large- argument expansion of the confluent hypergeometric function. 



4. As a nontrivial example, take the integral of the watermelon diagram treated in [Q only in D = 2 dimensions: 

D ^ fl\ ^ f d^p, d°p2 d^Pz 1 1 1 1 , . 

7 (2^)^(27r)^(27r)^l + pf l + p2l + pil + (pi+p2 + p3)2- ^ 

This integral has the powers 

fli = ^2 = ^3 = ^4 = 1, (19) 

and we identify the line momenta as 

qi^Pi, 92=P2, Qs^Ps, q4=Pi+P2+P3, (20) 

such that the matrix M is 

fli + 04 04 04 

M = 2 I 04 02 + 04 04 I , (21) 

04 04 03 + 04 

det M — 8(010203 + 010204 + 010304 + 02O3O4). (22) 
For the function I^^{k), we then obtain in any dimension D 
1 



^ ~ 24:D^3D/2 

^ 9D{2 + D)k'^ 9D{2 + D){4 + D)k^ 3D {2 + D) {IOA8 + 522D + 81D^) 
^ 32 128 4096 

91? (2 + D){A + D) (2576 + 9181? + lllD^) 
40960 

D{2 + D) (564864 + 397744D + 1109161)2 + I5228i:»3 + 891£|4) 
^ 65536 
3L> (2 + L>) (4 + D) (29651840 + 15696528i:' + 3452148L»2 + 3910681)3 + 19683i:'4) 

9175040 

31? (2 + L>) (4 + D) (1419854080 + 843338336L' + 2125088401)2 + 295623001)3 + 23449501)^^ + 857791)^) 



83886080 



For D = \ this reduces to 



(23) 
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27k^ 135^3 14859^4 97497k5 3268929^6 63271629k^ 22569248565k8 
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262144 



16777216 



+ 0(^9) 



and for D — 2 to 



28773 



9k2 27k3 453^4 1647^5 15157k^ 157293^^ 3720699kS 
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thus extending easily the expansions in Q| . 
For D = 3, the expansion reads 
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2127^9/2 



135^2 945k3 150435^4 1206387^^ 48595005^^ 1079675235k^ 



+- 



32 128 
432899207685^8 
16777216 
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8192 



65536 



262144 



(24) 



(25) 



(26) 



5. Having developed the tools for finding perturbation expansions of Feynman integrals, it remains to study the large- 
order behavior, and to find suitable methods for the resummation of the expansions with high accuracy. Together with 
the automatized generation of the Feynman diagrams of Rcf . , this will open the way for an "industrial production" 
of high- loop expansions for critical exponents. 
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